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Abstract. In this paper we present high spatial resolu- 
tion radio images of the puzzling binary system (3 Lyrae 
obtained with MERLIN at 5 GHz. We find a nebula sur- 
rounding the binary with a brightness temperature of 
(11 000 ± 700) K approximately 40 AU across. This defini- 
tively confirms the thermal origin of the radio emission, 
which is consistent with emission from the wind of the 
B6-8II component (mass loss of order of 10 -7 MQyr -1 ), 
ionized by the radiation field of the hotter companion. 
This nebula, surrounding the binary, is the proof that f3 
Lyrae evolved in a non-conservative way, i. e. not all the 
mass lost by the primary is accretted by the secondary, 
and present measurements indicate that almost 0.015 Mq 
had been lost from the system since the onset of the Roche 
lobe overflow phase. 

Moreover, the nebula is aligned with the jet-like struc- 
tures inferred from recent optical measurements, indicat- 
ing a possible connection among them. 
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ncnts is indicated by stationary optical and UV emission 
lines (Batten & Sahade 1973; Hack et al. 1975| ) as well as 
from the analysis of UV light curves (Kondo et al. 1994). 
Although (3 Lyrae has been observed extensively at other 
wavelengths for decades, radio observations are still too 
scarce to form a clear picture of its radio properties. The 
source was firstly detected by Wade & Hjellming ( 1972Q at 
2.7 and 8.1 GHz. During a successive monitoring program 
(Gibson 1975) the source always exhibited a spectrum con- 
sistent with thermal radio source with a radius ~ 50 AU, 
if a temperature of typical of an H II reg ion (10 4 K) is as- 
sumed. However, Wright & Barlow ( 1975 ) have noted that 
the observed slope of the radio spectrum (a=0.96) is in- 
termediate between that expected for a simple H II region 
and a stellar wind, implying that the physics underlying 
the radio emission of (3 Lyrae is more complicated than 
assumed in either of these models. 

In an attempt to distinguish the various components 
of (3 Lyrae which may be contributing to the radio emis- 
sion, we obtained high-spatial resolution radio maps of (3 
Lyrae using MERLIN. In this paper we describe the results 
of these observations and briefly discuss their implications. 



1. Introduction 



2. Observations and Results 



(3 Lyrae (3.4 < m v < 4.4, P = 12.91 d ) is a non- 
degenerate, semi-detached interacting binary syst em, lo - 
cated at a distance of 270±38pc (Perryman et al. 1997 ). 
The primary component is a B6-8 II star in contact with 
its Roche lobe resulting in rapid mass transfer (5 x 
10~ 5 M ( 7)yr~ 1 ) to its massive, unseen companion (Har- 



manec 1990 ). The nature of the secondary has been the 
subject of much speculation, from a flattened B5 star to 
a black hole with accretion disk. A B0V star embedded 
in a geometrically and optically thick accretion disk or 
a circumstellar shell, is currently the more favoured hy- 
pothesis (Hubeny & Plavec 1991, Harmanec 1992). The 



We observed (3 Lyrae at 6 cm (4.994 GHz) using the Multi- 
Element Radio Linked Interferometer Network (MER- 
LIN)Q array on 1996 December 8, 9 and 16, from 04:00 
to 22:00 UT. 

The observations were performed in phase referencing 
mode, which is the standard MERLIN observing mode for 
weak sources. This technique allows a better calibration of 
phase variation by observations of a bright phase calibra- 
tor interleaved with the observations of the target source. 
The phase calibrator 1846+322 was used and the flux den- 



presence of circumstellar plasma surrounding both compo- 



1 MERLIN is a national facility operated by the University 
of Manchester at the Nuffield Radio Astronomy Laboratories, 
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Fig. 1. The 6cm radio map oi (3 Lyrae. The contour plot 
shows the results of three separate 18 hour observations of 
(3 Lyrae combined into a single image. Levels of 0.18 (3a), 
0.24, 0.30, 0.36, 0.42, 0.48 and 0.54 mJy are shown. The 
position measured by Hipparcos is indicated with a large 
cross. The small inset ellipse is the half-power contour of 
the beam (60 x 47mas). 

Table 1. Parameters of the radio nebula around (3 Lyrae . 



Parameter 


Value 


Flux density 


3.2 ±0.2 mJy 


Size (FWHM) 


(145 ±12) x (100 ±8) mas 


Orientation (major axis) 


156.5° ±4° 


Radio position (J2000) 




RA 


18 h 50 m 04':800 ± 0:0060 


Dec 


±33°21'45'.'554±0'.'004 


Epoch 


1996.9 


Optical position" (J2000) 




RA 


18 h 50 m 04!794 ± 0!0002 


Dec 


±33°21'45'.'615±0'.'003 



a Hipparcos position reduced to same epoch 
including proper motion. 



sity scale was determined by observing 3C286, whose flux 
at 4.994 GHz is assumed to be 7.086 Jy. 

Data were edited and amplitude calibrated using the 
OLAF software package. The phase calibration, map- 
ping process and successive analysis were performed using 
the NRAO, Astronomical Image Processing System. To 
achieve the highest possible signal-to-noise ratio the map- 
ping process was performed using natural weighting and 
the CLEAN algorithm was used to extract information as 
close as possible to the theoretical noise limit. 

We detected the source at all three epochs, correspond- 



N-S direction. The noise level in the maps measured from 
several areas of blank sky is 6 x 10 _5 Jy/beam, which is 
consistent with the expected theoretical noise. There are 
no significant differences between the maps observed at 
each epoch and so we combined the data from the three 
observing runs to form the radio map shown in Fig. 1. 

The source morphology is dominated by a compact 
structure, aligned approximately N-S, probably (< 3cr 
level), embedded in a more diffuse nebula. The profile of 
the radio nebula along its principal axes is shown in Fig. 2. 
From the figure is evident that most of the 6 cm emission 
comes from the central region, with negligible contribu- 
tion from the more diffuse nebula, thus proving that a 
two dimensional Gaussian fit (dotted line) to the central 
part gives a fair description of the radio source. 

The position, flux density and angular size of the 
source are thus derived by fitting the two dimensional 
Gaussian brightness distribution to the map. The val- 
ues and their uncertainties, derived following Fomalont 
( |l989l ), are given in Table 1. The optical position of (3 
Lyrae reported by Hipparcos, allowing for the proper mo- 
tion, is also given. Differences between radio and optical 
positions (Aa = 0.075", AS — —0.061") may be probably 
related to the overall position error to be associated to 
the source used as phase calibrator, reported to be of the 
order pf 0.1" (Wilkinson et al. 1998). 

MERLIN is sensitive to structures up to 3.5" at 6 cm. 
No structure is seen outside of the nebula up to these 
scales, which would suggest that all of the flux from the 
nebula is detected in these maps. The difference between 
the flux measured here and previous measurements (e.g. 
Mutel et al. |l985| ; Leone, Trigilio & Umana |l994| ) is not 
surprising as (3 Lyrae is known to show variability at radio 
wavelengths. 



3. Discussion and Conclusions 

The radio observations here discussed are direct confirma- 
tion of a large plasma cloud surrounding (3 Lyrae due to 
mass loss from the system. The angular size of the neb- 
ula corresponds to a linear size of approximately 40 AU, 
much larger than the binary system itself (~ 0.5 AU). As 
the source has been resolved we can derive a brightness 
temperature Tb as follows: 



Tb ^1.36x10^ k = 11000± 700 K 

c'lC'2 



(1) 



where S is the flux in mJy at the 5 GHz, A is the observ- 
ing wavelength in cm and Q\ and 02 are the angular size 
of source in milliarcseconds. The derived brightness tem- 
perature implies a thermal origin for the radio emission. 

It can be argued that the source has, in reality, a much 
more complex morphology that the present MERLIN ob- 
servations can not resolve. If, for example, most of the flux 
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Fig. 2. Cross sections of the radio nebula along its prin- 
cipal axes (solid lines) together with the best fitting two 
dimensional Gaussian intensity distribution (dotted line) . 



model, we are assuming, will be no suitable. Still, in the 
former case, a higher brightness temperature would imply 
a non-thermal mechanism operating in the system. The 
presence of non thermal component has been inferred in 
about 30% of the early type radio emitting stellar systems 
(Bieging et al. 1989). All of these stars always exhibited 
spectral indices a < (S v oc v a ), while past (Gibson, 
1975) and more recent (Umana et al. 2000) radio obser- 
vations of (3 Lyrae point out a thermal spectrum [a > 0) 
for the system. 

Jameson & King (1978) modelled the radio source as- 
sociated to j3 Lyrae as a HII region. However, there is 
ample evidence, from both optical (Etzel & Meyer, 1983) 
and ultraviolet spectroscopy (Mazzali et al. 1992), of the 
presence of stellar wind in one or both components of the 
system, which can be responsible of the observed radio 
emission. 

Mazzali et al. ( 1992]) have developed a two- wind model 
to explain the peculiar UV spectrum of (3 Lyrae , consist- 
ing of " superionized" resonance lines, typical of hot-star 
winds, plus lines of moderately ionized species. The ob- 
served ultraviolet emission features can be modelled al- 
lowing the coexistence of two winds with different dynam- 
ical properties: a fast, but tenuous wind, associated to the 
hotter secondary component (Mh = 5 x 10~ 8 A/Qyr~ 1 , 
«h = 1470 km s _1 ) and a denser, but slower wind, as- 
sociated to the cooler B6-8II primary (M c = 7.7 x 




derive d by Panagia & Felli ( 1975 ) and Wright & Barlow 
( 1975 ). They showed that the observed radio flux (S) is 
related to the dynamical parameters of the wind through 
the relation: 



(2) 



where M and represent the mass-loss rate and the 
terminal velocity of the wind. 

We may ask under which conditions the thermal radio 
emission from the winds of the two early-type components 
of a binary system has the same radio properties as a single 
symmetric wind and thus the relation (j^) can be applied. 

The effects of binarity on thermal radio emission from 
early-type systems have been recently studied by Stevens 
(1995). Effects of possible interactions between winds, 
such as extra source of emission due to the shocked gas, 
are mostly function of the momentum ratio of the two 
winds and they would be negligible in systems, where one 
component has a dominant wind. The two-wind model 
outlined by Mazzali et al. (1992) foresees a case of a sys- 
tem when only one wind is dominant. Moreover, since the 
winds of the two stars will contribute to the measured flux 

according to the ratio of the relative factors (—), allow- 

f . . . °° . 

ing for difference in mass-loss and wind terminal velocity, 

this indicates that less than 0.5 % of the observed flux is 
attributable to the BOV star's wind. 

There is a class of stellar objects, the symbiotics, whose 
radio emission is interpreted in terms of a binary model, 
and where one of the two components has a dominant 
wind. Taylor & Seaquist ( |1984| ) attributed the observed 
radio emission to the stellar wind of the cool component 
of the system, which is ionized by the ultraviolet flux of 
the hot companion. They also showed that the geometry of 
the ionized radio emitting region of the intersystem mate- 
rial is determined by a single parameter, which is function 
of the physical parameters of the binary: separation (a), 
Lyman continuum luminosity of the hot component (L uv ), 
mass loss rate (M) and velocity (v) of the dominant wind. 
If the wind is completely ionized, the ionized region shows 
the same radio properties as a spherically symmetric ion- 
ized wind. In the following we will check if this scenario is 

suitable also for (3 Lyrae . 

By applying equation 14 from Taylor & Seaquist ( 1984| ) 
and assuming a binary separation of 61 solar radii (Har- 
manec, 1990), we obtain that the wind of the B6-8II is 
completely ionized if 



L uv > 5.8 x 10 62 ( — r 2 



1.8 x 10 45 photons s _1 (3) 



where the values of M and v for the B6-8II as deter- 
mined by Mazzali et al. (1992) have been adopted. This 
requirement is well satisfied by the flux of Lyman con- 
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contrary of the B6-8 II, and we can safely apply relation 

This formula has been derived for the case of spher- 
ically symmetric winds. However, Schmid-Burgk ( 1982 ) 
showed that the mass loss rate derived from radio flux 
densities is quite insensitive to source geometry and only 
in very extreme cases of deviation from spherical symme- 
try a geometry-dependent correction factor must be ap- 
plied. Thus even if the source morphology appears to be 
slightly elongated (axis ratio 1.45), we can use them quite 
confidently. 

We, therefore, can estimate the mass loss by using the 
relation: 

M = 6.7 x lO-SooS^Dj^ x g s )~ - 5 M^ 1 (4) 

where full ionization and cosmic abundances have been 
assumed, Voa is the terminal velocity of the wind, Se cm the 
observed flux density at 6 cm in mJy, Z?kpc the distance 
of the system, in kpc. gs represents the free-free Gaunt 
factor that, following Leitherer & Robert (1991), can be 
approximated with: 

T 3/2 

g s = 9.77(1 + 0.13 log— — ) (5) 

where T, in Kelvin, is the wind temperature. 

By assuming the stellar wind velocity of the B 6-8II 
component (v ~ 400 km s _1 , Mazzali et al. 1992| ) and 
constant temperature equal to the brightness tempera- 
ture, derived by the present observations (T = Tb), we 
obtain a mass loss rate of 



M = 5.6 x 10" 



M yr \ 



(6) 



Th is res ult is close to the value obtained by Mazzali 
et al. ( |1992| ) for the B6-8II component by modelling UV 
spectral lines, i. e. 7.16 x 10~ 7 A/Qyr -1 and indicates that 
the wind associated to the primary is predominant in the 
replenishment of the circumsystem material and consti- 
tutes the bulk of the system's mass loss. 

(3 Lyrae is approaching the end of a phase of strong 
mass transfer between components, which started 26000 
years ago (De Greve & Linnell 1994 ). The presence of such 
a nebula is the proof that the system had evolved in a non 
conservative way, i. e., with a fraction of the material from 
the loser lost by the system. This was already fo reseen 
by the evolutionary studies of De Greve & Linnell ( 1994 ) 
since their conservative solution lead to several discrepan- 
cies with the observed characteristics of the system. 

If we assume that during this phase of mass transfer 
the mass-loss rate was constant (5.6 x 10 _7 MQyr _1 ) we 
estimate that almost 0.015 Mq should have been lost by 
the system. 

High resolution optical interferometry combined with 



structures perpendicular to the orbital plane which are 
responsible for the bulk of the Ha and He I 6678 emiss ion. 
This scenario, already proposed by Harmanec ( 1992 ), is 
further supported by recent spectropolarimetric observa- 
tions of Hoffman et al. (199S), who suggested the presence 
of a bipolar outflow in (3 Lyrae , perpendicular to the or- 
bital plane (P.A. = 164°) and probably associated to the 
accretion disk around the mass-gaining component. 

The radio image presented in this work gives a hint of 
structure within the nebula, even if the obtained angular 
resolution did not allow us to trace the structure of the 
star's circumstellar matter at very small scale and thus 
to verify the association of the radio emitting features to 
the bipolar optical jets. One piece of information on this 
direction, provided by the present results, is that the radio 
nebula is almost exactly aligned with these optical jet-like 
structures (P.A. = 156.5° ±4°). 

It may be more accurate to picture the outflow as a 
stellar wind collimated by the thick accretion disk. The 
appearance of the circumbinary material would thus con- 
sist of a more extended component, due to the mass- loss 
from the system as a whole, plus an inner more compact 
component, related to the accretion disk created in the 
process of mass transfer from the B star to the hidden 
companion. 

To establish the feasibility of this scenario high resolu- 
tion radio observations, at frequencies higher than 5 GHz, 
aimed to probe the morphology of the inner structure of 
the radio nebula are necessary. This observations, by al- 
lowing to look more inside the binary's outflow, would pro- 
vide further insight into the formation of the circumstellar 
matter and arc of particular importance in understanding 
the connection of the observed radio nebula to the bipo- 
lar optical jets reported by Harmanec et al. (1996) and 
Hoffman et al. (|1998|). 
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